The Spectre of Too Many Species

Recent simulation studies examining the performance of Bayesian species delimitation as implemented in the BPP program have suggested that BPP may detect population splits but not species divergences and that it tends to over-split when data of many loci are analyzed. Here, we confirm these results and provide the mathematical justifications. We point out that the distinction between population and species splits made in the protracted speciation model (PSM) has no influence on the generation of gene trees and sequence data, which explains why no method can use such data to distinguish between population splits and speciation. We suggest that the PSM is unrealistic as its mechanism for assigning species status assumes instantaneous speciation, contradicting prevailing taxonomic practice. We confirm the suggestion, based on simulation, that in the case of speciation with gene flow, Bayesian model selection as implemented in BPP tends to detect population splits when the amount of data (the number of loci) increases. We discuss the use of a recently proposed empirical genealogical divergence index (gdi) for species delimitation and illustrate that parameter estimates produced by a full likelihood analysis as implemented in BPP provide much more reliable inference under the gdi than the approximate method PHRAPL. We distinguish between Bayesian model selection and parameter estimation and suggest that the model selection approach is useful for identifying sympatric cryptic species, while the parameter estimation approach may be used to implement empirical criteria for determining species status among allopatric populations

Digits Lost or Gained? Evidence for Pedal Evolution in the Dwarf Salamander Complex (Eurycea, Plethodontidae)

Change in digit number, particularly digit loss, has occurred repeatedly over the evolutionary history of tetrapods. Although digit loss has been documented among distantly related species of salamanders, it is relatively uncommon in this amphibian order. For example, reduction from five to four toes appears to have evolved just three times in the morphologically and ecologically diverse family Plethodontidae. Here we report a molecular phylogenetic analysis for one of these four-toed lineages – the Eurycea quadridigitata complex (dwarf salamanders) – emphasizing relationships to other species in the genus. A multilocus phylogeny reveals that dwarf salamanders are paraphyletic with respect to a complex of five-toed, paedomorphic Eurycea from the Edwards Plateau in Texas. We use this phylogeny to examine evolution of digit number within the dwarf−Edwards Plateau clade, testing contrasting hypotheses of digit loss (parallelism among dwarf salamanders) versus digit gain (re-evolution in the Edwards Plateau complex). Bayes factors analysis provides statistical support for a five-toed common ancestor at the dwarf-Edwards node, favoring, slightly, the parallelism hypothesis for digit loss. More importantly, our phylogenetic results pinpoint a rare event in the pedal evolution of plethodontid salamanders

The Probability of a Gene Tree Topology within a Phylogenetic Network with Applications to Hybridization Detection

Gene tree topologies have proven a powerful data source for various tasks, including species tree inference and species delimitation. Consequently, methods for computing probabilities of gene trees within species trees have been developed and widely used in probabilistic inference frameworks. All these methods assume an underlying multispecies coalescent model. However, when reticulate evolutionary events such as hybridization occur, these methods are inadequate, as they do not account for such events. Methods that account for both hybridization and deep coalescence in computing the probability of a gene tree topology currently exist for very limited cases. However, no such methods exist for general cases, owing primarily to the fact that it is currently unknown how to compute the probability of a gene tree topology within the branches of a phylogenetic network. Here we present a novel method for computing the probability of gene tree topologies on phylogenetic networks and demonstrate its application to the inference of hybridization in the presence of incomplete lineage sorting. We reanalyze a Saccharomyces species data set for which multiple analyses had converged on a species tree candidate. Using our method, though, we show that an evolutionary hypothesis involving hybridization in this group has better support than one of strict divergence. A similar reanalysis on a group of three Drosophila species shows that the data is consistent with hybridization. Further, using extensive simulation studies, we demonstrate the power of gene tree topologies at obtaining accurate estimates of branch lengths and hybridization probabilities of a given phylogenetic network. Finally, we discuss identifiability issues with detecting hybridization, particularly in cases that involve extinction or incomplete sampling of taxa

The role of historical and contemporary processes on phylogeographic structure and genetic diversity in the Northern Cardinal, Cardinalis cardinalis

Background Earth history events such as climate change are believed to have played a major role in shaping patterns of genetic structure and diversity in species. However, there is a lag between the time of historical events and the collection of present-day samples that are used to infer contemporary population structure. During this lag phase contemporary processes such as dispersal or non-random mating can erase or reinforce population differences generated by historical events. In this study we evaluate the role of both historical and contemporary processes on the phylogeography of a widespread North American songbird, the Northern Cardinal, Cardinalis cardinalis. Results Phylogenetic analysis revealed deep mtDNA structure with six lineages across the species\u27 range. Ecological niche models supported the same geographic breaks revealed by the mtDNA. A paleoecological niche model for the Last Glacial Maximum indicated that cardinals underwent a dramatic range reduction in eastern North America, whereas their ranges were more stable in México. In eastern North America cardinals expanded out of glacial refugia, but we found no signature of decreased genetic diversity in areas colonized after the Last Glacial Maximum. Present-day demographic data suggested that population growth across the expansion cline is positively correlated with latitude. We propose that there was no loss of genetic diversity in areas colonized after the Last Glacial Maximum because recent high-levels of gene flow across the region have homogenized genetic diversity in eastern North America. Conclusion We show that both deep historical events as well as demographic processes that occurred following these events are critical in shaping genetic pattern and diversity in C. cardinalis. The general implication of our results is that patterns of genetic diversity are best understood when information on species history, ecology, and demography are considered simultaneously

The role of historical and contemporary processes on phylogeographic structure and genetic diversity in the Northern Cardinal, Cardinalis cardinalis

<p>Abstract</p> <p>Background</p> <p>Earth history events such as climate change are believed to have played a major role in shaping patterns of genetic structure and diversity in species. However, there is a lag between the time of historical events and the collection of present-day samples that are used to infer contemporary population structure. During this lag phase contemporary processes such as dispersal or non-random mating can erase or reinforce population differences generated by historical events. In this study we evaluate the role of both historical and contemporary processes on the phylogeography of a widespread North American songbird, the Northern Cardinal, <it>Cardinalis cardinalis</it>.</p> <p>Results</p> <p>Phylogenetic analysis revealed deep mtDNA structure with six lineages across the species' range. Ecological niche models supported the same geographic breaks revealed by the mtDNA. A paleoecological niche model for the Last Glacial Maximum indicated that cardinals underwent a dramatic range reduction in eastern North America, whereas their ranges were more stable in México. In eastern North America cardinals expanded out of glacial refugia, but we found no signature of decreased genetic diversity in areas colonized after the Last Glacial Maximum. Present-day demographic data suggested that population growth across the expansion cline is positively correlated with latitude. We propose that there was no loss of genetic diversity in areas colonized after the Last Glacial Maximum because recent high-levels of gene flow across the region have homogenized genetic diversity in eastern North America.</p> <p>Conclusion</p> <p>We show that both deep historical events as well as demographic processes that occurred following these events are critical in shaping genetic pattern and diversity in <it>C. cardinalis</it>. The general implication of our results is that patterns of genetic diversity are best understood when information on species history, ecology, and demography are considered simultaneously.</p

An integrative approach to discovering cryptic species within the Bemisia tabaci whitefly species complex

Bemisia tabaci is a cryptic whitefly-species complex that includes some of the most damaging pests and plant-virus vectors of a diverse range of food and fibre crops worldwide. We combine experimental evidence of: (i) differences in reproductive compatibility, (ii) hybrid verification using a specific nuclear DNA marker and hybrid fertility confirmation and (iii) high-throughput sequencing-derived mitogenomes, to show that the “Mediterranean” (MED) B. tabaci comprises at least two distinct biological species; the globally invasive MED from the Mediterranean Basin and the “African silver-leafing” (ASL) from sub-Saharan Africa, which has no associated invasion records. We demonstrate that, contrary to its common name, the “ASL” does not induce squash silver-leafing symptoms and show that species delimitation based on the widely applied 3.5% partial mtCOI gene sequence divergence threshold produces discordant results, depending on the mtCOI region selected. Of the 292 published mtCOI sequences from MED/ASL groups, 158 (54%) are low quality and/or potential pseudogenes. We demonstrate fundamental deficiencies in delimiting cryptic B. tabaci species, based solely on partial sequences of a mitochondrial barcoding gene. We advocate an integrative approach to reveal the true species richness within cryptic species complexes, which is integral to the deployment of effective pest and disease management strategies

Molecular Systematics of the Deep-Sea Hydrothermal Vent Endemic Brachyuran Family Bythograeidae: A Comparison of Three Bayesian Species Tree Methods

Brachyuran crabs of the family Bythograeidae are endemic to deep-sea hydrothermal vents and represent one of the most successful groups of macroinvertebrates that have colonized this extreme environment. Occurring worldwide, the family includes six genera (Allograea, Austinograea, Bythograea, Cyanagraea, Gandalfus, and Segonzacia) and fourteen formally described species. To investigate their evolutionary relationships, we conducted Maximum Likelihood and Bayesian molecular phylogenetic analyses, based on DNA sequences from fragments of three mitochondrial genes (16S rDNA, Cytochrome oxidase I, and Cytochrome b) and three nuclear genes (28S rDNA, the sodium–potassium ATPase a-subunit ‘NaK’, and Histone H3A). We employed traditional concatenated (i.e., supermatrix) phylogenetic methods, as well as three recently developed Bayesian multilocus methods aimed at inferring species trees from potentially discordant gene trees. We found strong support for two main clades within Bythograeidae: one comprising the members of the genus Bythograea; and the other comprising the remaining genera. Relationships within each of these two clades were partially resolved. We compare our results with an earlier hypothesis on the phylogenetic relationships among bythograeid genera based on morphology. We also discuss the biogeography of the family in the light of our results. Our species tree analyses reveal differences in how each of the three methods weighs conflicting phylogenetic signal from different gene partitions and how limits on the number of outgroup taxa may affect the results

What are the consequences of combining nuclear and mitochondrial data for phylogenetic analysis? Lessons from Plethodon salamanders and 13 other vertebrate clades

<p>Abstract</p> <p>Background</p> <p>The use of mitochondrial DNA data in phylogenetics is controversial, yet studies that combine mitochondrial and nuclear DNA data (mtDNA and nucDNA) to estimate phylogeny are common, especially in vertebrates. Surprisingly, the consequences of combining these data types are largely unexplored, and many fundamental questions remain unaddressed in the literature. For example, how much do trees from mtDNA and nucDNA differ? How are topological conflicts between these data types typically resolved in the combined-data tree? What determines whether a node will be resolved in favor of mtDNA or nucDNA, and are there any generalities that can be made regarding resolution of mtDNA-nucDNA conflicts in combined-data trees? Here, we address these and related questions using new and published nucDNA and mtDNA data for <it>Plethodon </it>salamanders and published data from 13 other vertebrate clades (including fish, frogs, lizards, birds, turtles, and mammals).</p> <p>Results</p> <p>We find widespread discordance between trees from mtDNA and nucDNA (30-70% of nodes disagree per clade), but this discordance is typically not strongly supported. Despite often having larger numbers of variable characters, mtDNA data do not typically dominate combined-data analyses, and combined-data trees often share more nodes with trees from nucDNA alone. There is no relationship between the proportion of nodes shared between combined-data and mtDNA trees and relative numbers of variable characters or levels of homoplasy in the mtDNA and nucDNA data sets. Congruence between trees from mtDNA and nucDNA is higher on branches that are longer and deeper in the combined-data tree, but whether a conflicting node will be resolved in favor mtDNA or nucDNA is unrelated to branch length. Conflicts that are resolved in favor of nucDNA tend to occur at deeper nodes in the combined-data tree. In contrast to these overall trends, we find that <it>Plethodon </it>have an unusually large number of strongly supported conflicts between data types, which are generally resolved in favor of mtDNA in the combined-data tree (despite the large number of nuclear loci sampled).</p> <p>Conclusions</p> <p>Overall, our results from 14 vertebrate clades show that combined-data analyses are not necessarily dominated by the more variable mtDNA data sets. However, given cases like <it>Plethodon</it>, there is also the need for routine checking of incongruence between mtDNA and nucDNA data and its impacts on combined-data analyses.</p

Reticulate evolution: frequent introgressive hybridization among chinese hares (genus lepus) revealed by analyses of multiple mitochondrial and nuclear DNA loci

<p>Abstract</p> <p>Background</p> <p>Interspecific hybridization may lead to the introgression of genes and genomes across species barriers and contribute to a reticulate evolutionary pattern and thus taxonomic uncertainties. Since several previous studies have demonstrated that introgressive hybridization has occurred among some species within <it>Lepus</it>, therefore it is possible that introgressive hybridization events also occur among Chinese <it>Lepus </it>species and contribute to the current taxonomic confusion.</p> <p>Results</p> <p>Data from four mtDNA genes, from 116 individuals, and one nuclear gene, from 119 individuals, provides the first evidence of frequent introgression events via historical and recent interspecific hybridizations among six Chinese <it>Lepus </it>species. Remarkably, the mtDNA of <it>L. mandshuricus </it>was completely replaced by mtDNA from <it>L. timidus </it>and <it>L. sinensis</it>. Analysis of the nuclear DNA sequence revealed a high proportion of heterozygous genotypes containing alleles from two divergent clades and that several haplotypes were shared among species, suggesting repeated and recent introgression. Furthermore, results from the present analyses suggest that Chinese hares belong to eight species.</p> <p>Conclusion</p> <p>This study provides a framework for understanding the patterns of speciation and the taxonomy of this clade. The existence of morphological intermediates and atypical mitochondrial gene genealogies resulting from frequent hybridization events likely contribute to the current taxonomic confusion of Chinese hares. The present study also demonstrated that nuclear gene sequence could offer a powerful complementary data set with mtDNA in tracing a complete evolutionary history of recently diverged species.</p

The Application of DNA Barcodes for the Identification of Marine Crustaceans from the North Sea and Adjacent Regions

During the last years DNA barcoding has become a popular method of choice for molecular specimen identification. Here we present a comprehensive DNA barcode library of various crustacean taxa found in the North Sea, one of the most extensively studied marine regions of the world. Our data set includes 1,332 barcodes covering 205 species, including taxa of the Amphipoda, Copepoda, Decapoda, Isopoda, Thecostraca, and others. This dataset represents the most extensive DNA barcode library of the Crustacea in terms of species number to date. By using the Barcode of Life Data Systems (BOLD), unique BINs were identified for 198 (96.6%) of the analyzed species. Six species were characterized by two BINs (2.9%), and three BINs were found for the amphipod species Gammarus salinus Spooner, 1947 (0.4%). Intraspecific distances with values higher than 2.2% were revealed for 13 species (6.3%). Exceptionally high distances of up to 14.87% between two distinct but monophyletic clusters were found for the parasitic copepod Caligus elongatus Nordmann, 1832, supporting the results of previous studies that indicated the existence of an overlooked sea louse species. In contrast to these high distances, haplotype-sharing was observed for two decapod spider crab species, Macropodia parva Van Noort & Adema, 1985 and Macropodia rostrata (Linnaeus, 1761), underlining the need for a taxonomic revision of both species. Summarizing the results, our study confirms the application of DNA barcodes as highly effective identification system for the analyzed marine crustaceans of the North Sea and represents an important milestone for modern biodiversity assessment studies using barcode sequence